MOMENTUM REGULARITY AND STABILITY OF THE 
RELATIVISTIC VLASOV-MAXWELL-BOLTZMANN SYSTEM 



YAN GUO AND ROBERT M. STRAIN 



Abstract. In the study of solutions to the relativistic Boltzmann equation, 
their regularity with respect to the momentum variables has been an outstand- 
ing question, even local in time, due to the initially unexpected growth in the 
post-coUisional momentum variables which was discovered in 1991 by Glassey 
& Strauss |13j . We establish momentum regularity within energy spaces via a 
new splitting technique and interplay between the Glassey-Strauss frame and 
the center of mass frame of the relativistic collision operator. In a periodic box, 
these new momentum regularity estimates lead to a proof of global existence 
of classical solutions to the two-species relativistic Vlasov-Maxwell-Boltzmann 
system for charged particles near Maxwellian with hard ball interaction. 
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1. Introduction and formulation 

In 2003 it was shown for the first time that the full (Newtonian) Vlasov-Maxwell- 
Boltzmann system |17j has global in time unique classical solutions on the torus 
for initial conditions which are sufficiently close to the steady state. Then in 2006 
the rapid convergence on the torus for these solutions, as predicted by the 
H-theorem, was established. However it should be pointed out that this model 
is physically limited because it is not fully Lorentz invariant in the sense that 
the symmetries of the Maxwell system are inconsistent with those of the Newto- 
nian Boltzmann equation and serious difficulties are encountered in extending the 
method from [17] to the relativistic fully Lorentz invariant regime. In this work, as 
explained in the following developments, we overcome these difficulties and establish 
the momentum regularity as well as the existence of global in time classical solu- 
tions for the fully Lorentz invariant relativistic Vlasov-Maxwell-Boltzmann system 
with hard-ball interaction. 
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We study the following two-species relativistic Vlasov-Maxwell-Boltzmann sys- 
tem which describes the time evolution of charged particles: 

(LI) ' 



The collision operators, defined in (|1.4p below, are given by Q'^{F) = Q(F+,i^+), 
Q-(J^) = Q(F_,F_), Q±(F) Q(F+,^^_), and Qf{F) Q(F_,^^+). 
These equations are coupled with the Maxwell system: 



dtE -cV^x B = -47r 

dtB + cV^x E = 0, 
with constraints 



Jk3 1^ p\ pi J 



V, • = 4^ / {e+F+ - e-F^}dp, V, • B = 0. 

■/Ma 

The initial conditions are F±{Q,x,p) = Fo^±{x,p), E{Q,x) = Eo{x), and B{0,x) = 
Bq{x). Here F±{t,x,p) > are the spatially periodic number density functions 

for ions (+) and electrons (— ), at time t > 0, position x ~ {xi, X2, x^) £ T'^ == 
[— tTjTt]'^ and momentum p — (pi,p2,P3) S The constants ±e± and m± arc the 
magnitude of the particles charges and rest masses respectively. The energy of a 
particle is given by p^. = (m±c)^ -I- |pp and c is the speed of light. Note that 
here and below ± indicates two possible sign configurations. 

For number density functions i^-f (p) and -FL (p) a collision operator should satisfy 



Q{F+,F^){p)+ \ p \ Q{F^,F+){p)}dp^O. 





The same property holds for the other sign configurations. By integrating the 
relativistic Vlasov-Maxwell-Boltzmann system and plugging in this identity, we 
obtain the conservation of mass, total momentum and total energy for solutions as 



(1.2) 



dt 2 



( / p{m+F+{t) + m^F^{t)) + ^ [ E{t) x B{i)\ = 0, 

(m+p°i^+(i)+m_p"_F_(t)) + ^ / \E{t)\^ + \B{t)A^Q. 



The entropy of the relativistic Vlasov-Maxwell-Boltzmann system is defined as 

H{t) = - [ dxdp {F+it,x,p)\ogF+it,x,p) + F^it,x,p)\ogF^{t,x,p)} . 

Then the celebrated Boltzmann H-theorem for the relativistic Vlasov-Maxwell- 
Boltzmann system corresponds to the following formal statement 
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which says that the entropy of solutions is non-decreasing as time passes. 
The global relativistic Maxwellian (a.k.a. the Jiittner solution) is given by 

exp(-cp^/(/cBT±)) 



J±{p) = 



A'n:e±m\ckBT±K2{m±c^ / {kBT±)) ' 



where fcs > denotes Boltzmann's constant^ ^'2(') is the Bessel function K2{z) == 
e~^*(t^ — Vf'/'^dt^ and T± is the temperature. From the Maxwell system and 
the periodic boundary condition of E{t,x), we see that B{t,x)dx = 0. We 

thus have a constant B such that 

1 

|Tf3 



(1.3) / Bit,x)dx^B. 



Let [•, •] denote a column vector. We then have the following steady state solution 
to the relativistic Vlasov-Maxwell-Boltzmann system 

[F±it,x,p),Eit,x),B{t,x)] = [J±,0,S], 

which maximizes the entropy. 

We furthermore define the relativistic Boltzmann collision operator [5] as 

(1.4) Q{F^,G^)'^ f %l %l %W^\^[F^{p')G^{q')~F^{p)G^{q)]. 

It is written similarly for other sign configurations. Here the "transition rate" 
W^±|=F ~ W±\^{p,q\p' ,q') is defined as 

(1.5) W±\^ s a±|^(5, 9) S{pI + q^ - p'^ - q':^)S^'Hp + q - p' - ?')• 

The quantities here are s ~ s{p±,qz^), g = g{p±,q^)j ~ Sip±,q^) and <y±\zp{g,0) 
are defined exactly as in the following sub-section. 

The physical intuition provided by the Boltzmann H-theorem is to say that solu- 
tions should converge to their steady state, which is chosen by the initial conditions 
and the conservation laws (jl.2p . as time goes to infinity. Our goal in this work is to 
prove this global existence and rapid convergence in the context of perturbations. 

We define the standard perturbation f±(t,x,p) to J± as 

F±''=' J±+Vj±f±. 

We will plug this ansatz into (jl.ip to derive a perturbed system for f±{t,x,p), 
E{t,x) and B{t,x). The two relativistic Vlasov-Maxwell-Boltzmann equations for 
the perturbation / = [/+, /_] take the form 



with /(O, Xjp) = fo{x,p) ~ [/o,+(a;,p), /o.-(a;,p)]. The linear operator defined 
in (jl.Sp . and the nonlinear operator r-|-(/, /), defined in (|1.7p . are derived from an 
expansion of the Boltzmann collision operator (|1.4p . 

In particular, using (|1.5|) . we observe that the collision operator (|1.4p satisfies 

Q(J+, J+) = Q(J+, J_) - Q(J-, J+) - Q(J_, J_) = 0. 
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Then, with / = [/+, /_] and h = [h+, we can write the nonhnear operators as 

(1.7) r±(/, h) J±'^^Q{VJ±f±,VJ±h±) + J±'^^Q{VJ±f±,VJ^h^). 

Furthermore the hnearized cohision operators take the form 

(1.8) L±(/i) = -r±(/i,\/j)-r±(yj,/i). 

We estimate these operators in Sections [3] and H] 

In this hnearized system, the coupled Maxweh system is given by 

dtE ~cW,xB = -47tJ = -iTT f L+^Vj+.f+ - e_^\/j_/_l dp, 

(1.9) 7k3 I P+ P- J 

dtB + cV^x E = 0, 
with constraints 



(1.10) y^-E^inp^^Un [ \e+Vj+f+-e-Vj-f-}dp, V^-B^O 



The charge density and current density due to particles are denoted p and re- 
spectively. In computing p, we have used the normalization J±{p)dp = 

In all of the following developments, none of the physical constants 
vifill affect the results of our analysis. Then without loss of generality 
but for the sake of simplicity, we normalize all constants to one. We set 

J±{p)^J{p)^ 



An 



We drop the inessential ± notation all over, in particular we use the 
kernel notation a{g,0) — o'±|ip(.9,^) and we set p\. = p^ \/ ^ + \p\'^ ■> etc. 

Furthermore we assume that initially [Fo, Eq^Bq] has the same mass, total mo- 
mentum and total energy p.2p as the steady state [J, 0, B], then we can rewrite the 
conservation laws in terms of the perturbation [/, E, B] as follows: 

(1.11) / dxdpf+{t)Vj^ dxdp f.{t)Vj ^ 0, 

(1.12) / dxdpp{f+{t) + f^{t)}^7ipj^-^fdxE{t)xB{t), 

(1.13) / dxdpp'{f+{t) + f^{t)}Vj = -^fdx\E{t)\^ + \B{t)-B\'. 

We have used (|1.3p for the normalized energy conservation (jl.l3p . 

In the next sub-section we will discuss reductions of the collision operator (|1.4p . 

1.1. Collision operator in the Glassey-Strauss frame. In a pioneering work 
of Glassey and Strauss [14,, the collision operator Q was represented as follows: 

(1.14) Q{f,h)^ [ '-^i§2^B{p,q,Uj)[f{p')h{q')-f{p)h{q)]dujdq, 
where the kernel is 

(1.15) 5(p,g,^)- 
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Above the quantity s = s(p,q), which is the square of the energy in the "center of 
momentum" system, p + q ~ 0, is defined as 

The relative momentum, g = g{p, q), is denoted 



(1.16) g = v/2(pV-p-g-l). 

Notice that s ^ g^ + 4. Wc point out that this notation, which is used in [S], may 
differ from other authors notation by a constant factor. 
Tlie condition for elastic collisions is then given by 

^ p + q = p' + q'. 

In this expression, the post coUisional momentum are given as follows 

^-^ j^gN P' = P + a{p,q,uj)uj, 

q' = q~ a{p,q,uj)u}, 

where 

2b" + g>V{^-(^-^)} 



aip, q, w) = 



The Jacobian for the transformation {p,q) — > {p',q') in these variables [13] is 
^ ■ ' d{p,q) p^qO ■ 

Now we turn to a discussion of the collision kernel cr{g,6) in (|1.14[) . The kernel 
cr((7, 6) measure's the interactions between particles. See [SJin] for a physical discus- 
sion of general assumptions. We use the following hypothesis. 

Hypothesis on the collision kernel: We consider the "hard ball" condition 

cr((7, 6) = constant. 

This condition is implicitly used throughout the rest of the article. In fact to sim- 
plify the notation, without loss of generality, in the following we use the normalized 
condition a(g,6) = 1. The Newtonian limit, as c — > oo, in this situation is the 
Newtonian hard-sphere Boltzmann collision operator |24j . 

In the next section we will discuss our main results. 

2. Statement of the main results 

Let the multi- indices 7 and /3 be 7 = [7", 7^, 7^, 7^^] and /3 = [/3-^ , , /J"^] . Wc will 
use the following notation for a high order derivative 

— '^Xi ^X2 "X3 ^pi ^P2 ■ 

Wc sometimes also use the notation /3o, /3i, /32 and a to denote multi- indices with 
three components such as /3. Then if each component of /3 is not greater than that 
of /3i, we denote this by /3 < /3i, also /3 < /3i means /? < /3i and \f3\ < |/3i| where 
\/3\ ~ (3^ + (3"^ + (3^ as usual. We also denote a collection of weight functions by 
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Given a solution [f±{t,x,p), E{t,x), B{t,x)] to the relativistic Vlasov-Maxwell- 
Boltzmann system (|1.6|) . (|1.9p . and (|1.10[) . we define the full instant energy func- 
tional to be a continuous function, £N^e{t), which satisfies: 

(2.1) 8NAt)^ E \\md;f{t)f+ E mE{t),B{t)W. 

h\ + \l3\<N h\<N 

At time t = the time derivatives in SN.iiO) are defined customarily through 
equations (|1.6|) and (|1.9|) . In ()2.ip and the rest of this paper, the norm || • || denotes 
either the L'^if^ ^ ^p) norm or just the L^(T^) without ambiguity depending upon 
the variables in the functions being measured. Throughout the rest of this paper 
we furthermore assume > 4 and ^ > 0. 

We are now ready to state our main results: 

Theorem 1. Suppose that [/o,±, ^'Oi ^o] satisfies the constraint p.lOp and the the 
conservation laws (jl.lip . (|1.12p . (|1.13p initially. Fix £ > and N > A. Consider 
Fo^±{x,p) = J-t + \/j±fo,±{x,p)- There is a constant M > such that if 

SnAO) < M, 

then there exists a unique global solution [f± [t, x,p), E(t, x), B{t, x)] to the perturbed 
relativistic Vlasov- Maxwell- Boltzmann system (jl.6p . (jl.9p with (|1.10p satisfying 

Moreover F±{t,x,p) = J± + V^±f±{t,x,p) solves the relativistic Vlasov- Maxwell- 
Boltzmann system (jl.ll) . and F±(t,x,p) > if this is so initially. 

If £N+k.eiO) is sufficiently small for some fc > 0, then we have rapid decay as 

E \\w,d}f{t)f+ E W[E{t),B{t) - B]f <^N+kA^){l+t)-\ 

\l\ + \l3\<N \-i\<N 

where the constant B is defined in (|1.3p . 

There have been many investigations about various kinetic models for describing 
charged particles. Standard references include [3ll5llTl]. We refer to several results 
such as [UdlinKTnillllllHlllI]. In [TB], global classical solutions were constructed 
for the Vlasov-Poisson-Boltzmann system (no magnetic fields) via introduction of 
a nonlinear energy method for which the linear collision operator L is positive 
for solutions near Maxwellians. In [17j . an improvement of such a method led to 
construction of global solutions to the Vlasov-Maxwell-Boltzmann system in the 
presence of magnetic field. In j23j , such a construction was carried out in the whole 
space case using some new dissipation estimates. Even though the Vlasov-Maxwell- 
Boltzmann system can be viewed as a 'master system' for kinetic models, from 
general physical principle, the classical (non-relativistic) Boltzmann is not compat- 
ible with the (Lorentz invariant) Maxwell system, which obeys special relativity. 
It is therefore important to study the relativistic effects for the relativistic Vlasov- 
Maxwell-Boltzmann system (|l.ip and to generalize the result in [T7] to a relativistic 
setting. However, such a project was easily stopped due to a severe difficulty of lack 
of regularity in the momentum p variables for the relativistic Boltzmann equation. 

In the Glassey-Strauss frame (|1.14p and p.lSp . the following pointwise estimates 
were discovered by Glassey & Strauss [13] in 1991: 

(2.2) \Vj,q[\ + \Vj,p[\ <C{qf {l + \p- ^r/'l{|p..|>|px<.|3/2}) . 
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This is a sharp estimate at least in regards to the p growth. Furthermore 

\y,p',\ + \V,ql\<C{qfp°. 

Although in this second estimate, no attempt was made to further refine it. Notice 
that higher derivatives do not exhibit momentum growth in p: 

The q growth above does not cause any diSiculty because we have strong exponential 
decay in that variable in the linearized collision operator T. But the p momentum 
growth in ()2.2p . introduces high order growth of p in {|Vp(7,-| + |Vpp-|}^, within 
the highest order derivatives of dpf{p') with |/3| = N. 

Such a growth phenomenon is purely a relativistic effect, which is absent in the 
classical setting. Consequently, no regularity for the momentum variables, even 
local in time, for the solutions of the Boltzmann equation has been obtained. This 
is because of the presence of highest polynomial growth rate in p in the highest 
order derivatives of the solutions, which makes it impossible to close the estimates 
in any reasonable energy norm. Up to now, all existing work for the relativistic 
Boltzmann equation only involves spatial and temporal regularity. Unfortunately, 
it is necessary to obtain momentum regularity to study the relativistic Vlasov- 
Maxwell-Boltzmann system due to the Lorentz force term: {E + p x B) ■ Vpf. 

2.1. Collision operator in the center of mass frame. Our key idea to overcome 
such a severe difficulty in the Glassey-Strauss frame is to resort to the following 
center of mass representation of the Boltzmann collision operator. We remark that 
the study of such a center of mass frame was initiated recently in the absence of 
momentum derivatives For a function G : x x R3 ^ j^a ^ jg. 

sB{p,q,uj) ' '\ f .1 m " "^ 
(To G{p,q,p,q)= duj V0 G{p,q,p ,q 

where B{p, q, uj) is given by (|1.15|) and (p', q') on the l.h.s. are given by (|1.18|) . On 
the r.h.s. we use that = V0{p,q) is the M0ller velocity given by 

2-3 ^ o-Q-Q- 

2 p"q'^ 

The post-collisional momentum, on the r.h.s. can be written: 

„ dcf P + q Q f , . _,w , {p + q)-u 

„ dci p + q g f , , Ap + q)-u 



where p = [p^ + q^)/ ^/s. See [261 Corollary 5] for basic properties of such a center 
of mass frame. In particular the Jacobian (|1.19p effectively also works here as 

(2.5) / dpi dq f duJV,aig,e)Gip,q,p",q") 

JR3 JK3 JS2 

dp dq duj V0 a{g,6) G{p" , q" ,p, q). 

A more detailed explanation is given in [26l Corollary 5 and (23)]. 

Clearly, there is also a problem in the center of mass variables from (|2.4|) . In 
these variables it is straightforward to compute that high momentum derivatives 
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of p" and q" create high singularities when p = q and p = —q. These two distinct 
problems in each separate representation formula for the relativistic Boltzmann 
collision operator illustrate the main reason why it has remained an open problem 
to prove energy estimates with momentum regularity. 

To resolve these difRculties, and to prove the main nonlinear estimate in Theorem 
[2] below, we will split the desired estimate into two different cases. These cases 
correspond to the following two different integration regions: 

(2.6) ^'=^'{H<l}U{|p|>l,H^ <2g°}, A/='{H>l,Hi/™>2qO}. 

Here m » 1 is taken to be a given large positive integer. On the set A, we can use 
the Glassey- Strauss frame (|1.18p . Large growing polynomial momentum weights 
in p, as described above and in (j2.2|) can be controlled by the factor J^/*(g) since 
\p\ ^ {q^)^ on A. On the other hand, on the region A^, we will use the center- 
of-momcntum variables \2A\ . Note \p\ > 2\q\ so that |p ± q| > -^y- > i. Then the 
deficiency in these variables, namely that derivatives of (|2.4p create singularities 
(even though there is no momentum growth at infinity), is fortunately avoided on 
the region Ac] meaning that our estimates in this region are safe as well. Without 
such a magical use of the center of mass frame, it is still an outstanding open 
question if one can control the growth in p solely within the Glassey-Strauss frame. 

We would like to say that we think it would be interesting to study this analogous 
problem in the whole space E.^ including the convergence rates, as in [6j[7,i23,i27j . 

2.2. Notation. In addition to the notation from (|2.ip . we will use the spaces 

\M2= jf dx f dp\h{x,pW, \h\2 = jf dp\h{p)\^. 

Similarly any norm represented by one set of lines instead of two only takes into 
account the momentum variables. We also define \\h\\i, = \\h\\2, which is justified 
by (|4.ip later on. The L^(Rp) inner product is denoted (•, •). We use (•, •) to denote 
the L^i^l X Kp) inner product. Now, for ^ <S M, we consider the weighted spaces 

||/i||2,£ ||wf/i||2, \\h\\uj=^ \\wih\\„, \h\2,i \wih\2, \h\^^i'^= \weh\^. 

We will furthermore use A < B to mean that 3C > such that A < CB holds 
uniformly over the range of parameters which are present in the inequality (and 
that the precise magnitude of the constant is unimportant). The notation B > A 
is equivalent to A < B, and A ^ B means that both A < B and B < A. We 
additionally use C > to denote a generic positive large constant and c > to 
denote a small constant; their exact values are considered to be inessential. 

2.3. Organization of the paper. In Section[3]we will prove the estimates for the 
momentum derivatives of the non-linear collision operator (|1.7|) using the splitting 
into the Glassey-Strauss frame and the center of momentum frame. Then in Section 
|3]we will use these nonlinear estimates to deduce quickly several linear estimates, 
using also Lastly in Section [5j we show how to use our estimates to prove the 
global existence and rapid decay, following the arguments from [TU[5D]. 
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3. Momentum derivatives of the nonlinear collision operator 



Recalling the decomposition of A and Ac in ()2.6|) . consider the smooth test 
function x G ^'^"([O, oo)) such that < x < Ij and x{p) = 1 for P G [0, 1] with 
= for p > 2. We use the splitting 1 = Xa{p, q) + XA^ {Pi q) with 



XA{p,q)'^x{p") + {l-x{p"))x 



XA.(P,<Z) = (1-X(/)) 1-x 



We split r(/i,/2) = rA + r^, as 

sB{p,q,ij) 



(3.1) 



dujdq 



[h{p')f2{q') - h{p)f2{q)]XA{p,q), 



Here without loss of generality, we have taken /i and /2 to be scalar functions. 
Using these important decompositions, we will prove the main estimate: 

Theorem 2. We have the following nonlinear estimate for any > 0; 

\{wldpT{f^,f^),dpf:,)\<\dpf^\2,i \dpJi\2AdpJ2\2,i. 

Here we can include any i>Q. Then for I7I + < N with N > 4 we have 
{^^^;T{f^J,),^;f^)\< Wdpshj n E Il^ft/,ll2,.. 

i=l,2|7i| + |;3i|<Af 

The second estimate follows easily from the first and Sobolev embeddings. 

Now this theorem will follow directly from our Lemmas [2] and [5] below. It will 
be our focus in the rest of this section to prove these estimates. 

3.1. Estimates in the Glassey-Strauss Frame Fa- To avoid taking derivatives 
for the singular factor of |cj • — ^)| inside B{p,q,ijj) for dpTA in (|3.ip . we 
introduce the following change of variables q ^ u (for fixed p) as: 

(3.2) u = p°q- q°p. 

By p.2p . we have that q ^ + ^ and taking norms on both sides yields 



q° ^{u-p) + \/{u-pY + |u|2 + (p0)2_ 
Such a transformation p.2[) therefore defines an invertible mapping with 



dui 



du 



Since Ic.- (5^-^)1 



= det 



qjPt 



q' 







duj 
dqj 



(«,i = 1,2,3), 
{p'f 



q^ - ' - " q 
, we can express V a from as 



p-q) > 



(p°) 



0\2 



p"q 







(3.3) Fa = / dudu 



dq 



du 



sB \uj ■ u\ 
p°q" 



VJiq){fiip)f2{q) - fl{p)f2{q)}XA{p,q), 
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where now 



We take a high order derivative dp of p.3p to obtain 



< 



E 



|P|~<9' 



■E 



where the sum is over /3o + /^i + /32 < Furthermore 



Jiodu Kl |aft/i(p)9/5j2(g)|, 



(3.4) Kl=Ki{u,p,u:)''^\u:-u\ 



dq 



du 



sB 

^ J'^Hq)xAip,q) 



J-y\q). 



to 



Also yU^^ = w) is the term which results from applying the chain rule 

the post-coUisional velocities p' and q' . Here //^^ contains the sum of products of 
high order momentum derivatives of the smooth functions p' and q' . The next step 
is to reverse this change of variables p.2p to go from u back to an integration over 
q. After that change of variables: 

Lemma 1. On the set A, we have the following estimates 



< 



wl(p) \K^^,{p°q~q%,p,io) 
Similarly, we also have the upper bound of 

/^|(A-9Vp,^)|<(g>". 

Above n > I is a fixed large integer which depends upon £ >0, l3, /3q, and (32- 

Proof of Lemma\^ Wc start with the estimate for Clearly, up to constants, 
^^p\ is a sum of products of terms of the following form 

{dp,py-{dMr\ 

where (32, 71 and 72 are suitable multi- indices which are all < \(3\. (Note /3i and 
P2 in the previous display need not be the same as those in It is therefore 

sufhcicnt to estimate the size of these derivatives from above. This follows from 
the multi-dimensional generalization [4], from 1996, of the Faa di Bruno formula 
(1855). Consider the case |/3| ~ 1. After the change of variables p.2p we have 



P 



p + a{p, u,uj)uj, 
-n+-oP-aip,u,uj)uj, 



where 



a(p, u, uj) 



2{p° + q°) UJ -u 



{pO + qOy-{u;.{p + q)} 



_ N 
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Our goal will be to estimate derivatives of these functions. We thus compute 



— — = Okj + UJk- 

dpj 



dpj 



dp, (pO) 

We compute the final derivative as 

da _ I dN 
dpj D dpj 
dN 
dpj 

dp, 



pkPj 



' Pk , da{p,u,uj) 



dpj p° 



dpj 



2uj ■ u 



N dD 
^ D^d^' 

Pj_^d(t 
P^ dpj 
2{u ■ p)uj + 2pj 



^(M-p)2 + |u|2 + (pO)2 

The derivative of the denominator is further given by 

^0^ 



dpj 



dpj 

3 

2w ■ (p + g) ^ Wfc 



k=l 



rSkj 



Pk 



_d_ 

dpj \p' 



UkPj 



In particular we can write the whole derivatives of {p' , q') as 

dp'k ^ .fkj{p,q,(^,P°,q°,u) 

dpj 



dpj 



{p°)''Db{^{u ■ pY + |m|2 + [p^YY 
gkj{p,q,i^,p°,q",u) 



[p^YD^^Jiu-pY + W + {p'^YY ' 



where fkj and gkj are smooth polynomials in the variables {p,q,uj,p^,q'~',u). Fur- 
thermore, a, b and c are positive exponents which depend upon k, j and the form of 
(p', q'). These polynomials and exponents are quite lengthy to compute. However, 
the key observation in these calculations is that the denominators, e.g. p", D, and 
y/ {u ■ pY + |wp + {p'^Y ^-re in all cases uniformly bounded from below, so that no 
singularities are present. 

Therefore, after applying the reverse change of variables w — s- p^g — q'^p to p.2p 
we can see that we always have the crude upper bound (for some n > 0) of 



dp'k 


+ 


dq'k 


dpj 




dpj 



< 



Hp) {q)y 



This estimate will conclude the second estimate in Lemma [T] if all of the derivatives 
arc first order derivatives. Note that the exact value of n is in fact unimportant 
to our argument. The crucial observation now is that this pattern repeats for the 
higher derivatives of order /3 with > 1. In particular, we sec that 

fl3]{p,q,(^,P°,q°,u) 



dpp'k 

dpq'k 



(pOYDb{y/{u-pY + |up + (p^YY 

9i3j{p,q,(^,p°,q",u) 

ip"YD''iViu-pY + |wP + ip")'^Y' 
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where again fpj and gpj are smooth polynomials in the variables {p,q,uj,p'^,q'-' ,u). 
And once again a, b and c are (different) positive exponents which will depend 
upon /3, k and the form of {p' , q'). This form of these high order derivatives 
is quickly deduced from the standard rules of differentiation, and for instance a 
simple induction procedure. However to compute the exact expressions of fpj and 
gpj seems to be quite difficult. A key observation is that computing these high 
order polynomial expressions explicitly is in fact not-necessary to our argument. 

Again the crucial point is that the denominators of dpp'f^ and d/^q'/. are uniformly 
bounded from below by a positive constant, so that no singularities are present. 
This implies that, after applying the change of variables u — > p^q — q'-'p, the higher 
order derivatives are all similarly bounded above as 

\dpp'k\ + \df^q'k\<{{p) {q)y\ 
for some (different) n > 0. After that, the second estimate for /i^J in Lemma [T] is 
a consequence of the region A from (|2.6p . Wc conclude the estimate for y^t^^ . 

The estimate for ()3.4p is directly similar. The point is again that momentum 
derivatives in p.4p end up creating ratios of polynomials in the variables w, u, 
p, p", and q^. The exact expressions created by the high order derivative d/^g in 
p.4p is apparently quite difficult to compute in general, but fortunately this is not 
necessary. Instead it is easily seen that the denominators of these rational functions 
created by the high order derivative, dp^^ are all uniformly bounded from below. 
Thus after reversing the change of variables in (|3.2I) the term (|3.4[) must be bounded 
from above by a constant multiple of {{p) (g))" for some n > 0. □ 

Lemma 2. We have the following estimate for Ta with any £ > 0: 
\wUp)dprA\ < e-^l^l''" / dqdu; 1 x „ ji{q) 

(3.5) X {l(5/3i/l)(p')(5/../2)(9')l + l(5ft/l)(p)(a^./2)(9)|}. 

/3i+/32</3 

Moreover, from that estimate one can deduce the following uniform bound 

(3.6) \{w]dpr^if,j2),dpf3)\<\dpfs\2 J2 

Proof. The proof of p.Sp follows directly from the previous Lemma [1] by noting 

that J-^{q) > e'^IPl'^" > (1 + \p\ + |g|)" for any large n > since \p\ < (g^)™. 
To establish p.6p . we apply the Cauchy-Schwartz inequality to obtain 



\{w'jdi3TA,df3f3)\ 



< L / —ol^{\dpJlip')dMq')\ + \df,Jl{p)dpj2{q)\}\d^f3{p)\ 

J"'{q) 



<\dpf3\2 E j dujdqdp \{dpJ,){p'){dfU2W)\ 



2 



01+P2<I3 



+ I9/3/3I2 V \dl3jl\2\dl3^f2\2- 



I3i+I32<l3 



We finally complete the proof by making pre-post coUisional change of variables 
(p,g) ^ (y,q') from (inil). □ 
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3.2. Center of Momentum Frame. In this section we prove estimates for the 
term Ta^ from (j3.ip . We take /3 momentum derivatives of to obtain 



|p|l/™>q' 



o\dp,V,\^)\idpJl){p"){dpj2){q")4l\ 



Here Kq^ is the coUection of sums of products of momentum derivatives of p" and 



q", from (|2.4p . which result from the chain rule of differentiation. Again the sum 
is over the multi-indices Po + (3i + (32 < (3. We then have 

Lemma 3. Let |p| ™ > q^ with m large, as in (|2.6p . For some integer n > 1, which 
depends upon j3 ^ 0, we have the following estimates 

— — + \dpP I + Idpq \<{q) . 
To prove Lemma [3] we will use the following: 
Lemma 4. Let |p| m > with m large. Then for any /3 7^ 0, we have 



\dp9\ < 



1^1 



9 

{q)\^\ 



Ob - 



< 



9 



1 



< 



9^ 



Proof of Lemma^ We shall use an induction over |/3|. If \/3\ = 1, we have 

<M 



dp. I - 



\9p,9 

1 

9 

\dp,Vs\ 
1 



^ ( 

- —q - 

9 VP 



dp,9 



9' 
9dp,9 



< 



9dp,9 



< 



< 



9 



9 



l'9p.g| < M 
5^ '"5^ 



(g> 



\dp,9\ 



< 



Here we recall that s = g'^ + 4. Assume that the estimates are valid for /?. Now for 
(|/3| + 1)— th order derivatives, from the induction hypothesis for -, wc have 



\9f)dp,9\ 



/3i</3 



Ob 



d. 



Pi 

—q -qi 



^ 0#/3i</3 ^ 



< 



{q) 



1/31 + 1 



Here (and below) we use that g > 1 on \p\ >" > q'^ . This follows from 
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which itself is a consequence of the inequahty ^ 3 ^ b ^ 91 (see [111 Lemma 



p"q" 



3.1]) on the region \p\ ™ > (7° from p.6p . For the next step, similarly note that 



< 



< 



< 



/3i+/32</J 

E 

/3i+/32</3 



(g)l/3|-|/3i|-|/32| (g)l/32| + l 



This last estimate in particular holds because we have, by the induction assumption 
and a direct calculation, for any multi-index < a < /? that 



da 



< 



(9) 



Here again we used that 5 > 1 on \p\ >" > (?°. 

Simialrly for any multi-index satisfying < a < /3 we have that 



da 



1 



< 



\a\ 



\dag\ < max<^ 5 



(9> 



With that, we again use the induction hypothesis to obtain 



< 



< 



< 



da 



E 

ft+/32</3 



E ^^^^ i 5' ■ 

/3i+/32</3 
(g)l/3| + l 



(g)l/3|-|ft|-|/32| ^ (^qyi3i\ (^q)\l32\+l 



9 



For the last case, we do a similar calculation as 



gdpig 



1 



< 



/9i+/32+&+/34</3 

(q)lftl (5)1/321 (q)lfel 



E 



/3i-l-ft-l-/33+/34</3 



5 



5 



■ max < g 



(g)l/34| ^ (q)l,9| + l-|ft|-|/32|-|/33|-|/34| 



< 



(9) 



1/31 + 1 



In summary, the desired estimates follow by via the induction hypothesis. 
We now use the estimates from Lemma U] to prove Lemma [31 



□ 
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Proof of Lemma\3i We first show the decay of di3V0. Consider g, and recall p.7p 
Next recall = from (gS]). Then using Lemma |4] we have 



\dfiv, 



dp 



pOqO j 



< max < g 



I3i+P2<l3 



9 



9 



dpi9di32Vsdi3-p,-i3^ 

{p) 



pOqO 



< 



g2 pOqO 



This completes the estimate for a high-order derivative of V0. 

To show \di3p"\ + \df}q"\< (g)" , we note from jM]) for > that 



9fi9 



E 

a<l3 



, ( P + g) ■ 



We estimate each of these terms individually. The first term is trivially bounded. By 
LemmaSl for the second term on the right side of dfip'l we use \dpg\ < — < (g)''^' 
IpI™ Foi' the third and last term on the right side of dfip'l wc notice 



da {Pi + qi) 



(p + q) -uj 
\p + q\^ 



< 



\p + q 



The first inequality holds generally; the second inequality holds on \p\ ™ > g° from 
(|^ . Now for the term §(p - 1) we notice that |f (p- 1)| < P° + q° < (p) on 
bl" ^ The first inequality in the previous chain holds because of ^ < 1 and 
the definition of p from (|2.4p . We conclude that if /3 — a then the third term 
on the right side of djsp'- is bounded as in Lemma [3] since in this case |a| = > 0. 
It remains to estimate the last term on the right side of when |/3 — a| > 0. 

To this end, notice that |(p - 1) | ( ^^^^ - l) • Therefore by Lemma H 



a^(f(p-i) 



<\{dpg){p-l)\+ Yl 

1 



< 



l/3| 



E 



max 



E 1 5 

/3i#/3,/3i+/32</3 I 



9/3- p 1-132?° 



(q) 



(?) 



< 



(q) 



1/31+2 



This use several previous estimates, and holds for a general multi-index /3 7^ 0. 

Collecting all of these completes the estimate for 9/3pf in Lemma |3l Notice that 
the estimate for |9/3g-'| is exactly the same. □ 
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Lemma 5. Fix |/3| > 0. Then we have the uniform estimate 

(3.8) Id.TAj ^ dqdu; lb|V,.>,o I(gft/i)(y0(gfe/2)(g^0l 

+ E|9/Ji./i(p)| f dqdw l|p|i/^>,o J1mi,S} l^fe/2(g)l- 

^460116 i/ie sitm is over multi-indices /3o + /3i + /32 < /3. 
Moreover, for i > 0, we obtain 

(3.9) \{w'edprAAhJ2),d0fs)\<\dpf3\2j \dpjike\dpj2\2,i. 

Proof. First p.8p follows directly from the previous Lemma |3] and the fact that 
J~^/'*(q) > (g)" for any n > 0. We then have the upper bound of 



wldpT^^,dpf:,)\<Y, j wl{p)v,j''\q) \{dpJ,){p"){dpj2){q")dpfM\ 
+ f dqdpdu: wj{p)v,J^/\q) \dpj2{q)\ \dpj,{p)dpf3{p)\. 



Above the sum is over multi-indices /3i + /32 < The second term above clearly 
has the desired upper bound in p.9p using Cauchy-Schwartz. 

For the first "gain term", notice from [HI Lemma 2.2] that we have the estimate 
wtip) < we{p")we{q"). We remark that the estimate [HI Lemma 2.2] is true for 
any variables satisfying the conservation laws (|1.17l) . Now using Cauchy-Schwartz 
we obtain the upper bound (using also V0 < 4) 

\dpf3\2,e(^J2 J wUp''WM')\df3jiip'')df3j2iq'')fdpdqdu;^ ' 

< \di5f'i\2,i \ Y I V0 wj{p")w'j{q") \dpji{p")dp^f2{q")f dp'dqdu 
Once again the pre-post change of variables from (|2.5p establishes Lemma [5l □ 

4. The linear estimates 

We will use the estimates proven in the previous section for F from (|1.7p to prove 
the linear estimates in this section. Recalling (|1.8p and (|1.7p we write 

L{h) = [L+{h), L^{h)], L{h) = v{p)h - K(h). 

Here we recaU that /i = hJ\. From (|1.4p . (|1.7p and (|1.14p we can define 

r^""(/i,/2) = / f^^%^yj(^/i(y)/2(g'), 

r'°^^(/l,/2) = / d^dq '2E^PlSl^^)f,[p)f,[q). 



Above /i and /2 are scalar functions. Then following (|1.7p we can write 
(4.1) .(P) = 2r'-^ (1, x/j) = 2 / dcd, Jiq). 
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With these developments the operator K{h) = [K+{h), K-{h)] can be expressed as 

(4.2) i^±(/i) = rs""(/i±,yj) + rf""(/i±,VJ) + r± . 

Then from [25j Lemma 3.1] we clearly have that i'{p) « for Co- > 0. Fm-thermore 
Lemma 6. Let \P\ > 0, then \dpiy{p)\ < C {p)~'^ . 

Proof. As in (|4.ip . we apply p.Sp and ()3.8p with /i = 1 in the loss terms. □ 



Proposition 7. Lef |/3| > and fix £ > 0. For any small rj > 0, there exists a 
large R = R{ri) > and C ~ C{r]) > such that 

{w^,dpMp)h},dph)>\dph\l,^7^ - Cr,\l<Rh\l 

l"l<l/3| 

Here 1<r{p) is the indicator function of the ball of radius R centered at the origin. 

Proof. We will prove the desired coercivity estimate for a real valued function h 
to simplify notation; the result follows trivially for a vector valued function h = 
We expand out the inner product as 

{wjdp{iy{p)h},di3h) = / + Cf / wl dp^v dp-p^h dph 

= \drih\lj + XI '^ft / '^P "'^(^^) dp^v{p) dp-p,h{p) dph{p). 

0</3i</3 •^^'^ 

Here C^^ is the constant which results from the high order differentiation. Since 

\(3i\ > we have from Lemma |6] that \di3^v{p)\ < C {p)~^ . 
Then, for fixed i? > 0, we split the second term above as 

V Cf / dp wj dfi^v{p) df3-f3^h dph^ j + 

0</3i</3 -^KS "'|p|<-R J\p\>R 

On the unbounded part we use Cauchy-Schwartz as follows 

On the bounded region we use the compact interpolation of Sobolev-spaces 

/ d.p \dp^p,hdph\< / S E dp \dp^p,h\'^ + fi\dph\" 

J\p\<R 0<pi<l3 J\p\<R [o</3i</3 

< 77' V / dp + C^- f dp \h\^ 

i_i ,o,Jm<r J\p\<r 



<V J2 \d^h\li + Crj\l<Rh\ 



Above the different > variables are allowed to be arbitrarily small. This com- 
pletes the desired estimate. □ 
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Proposition 8. Let \/3\ > and £ > 0. For any small ?; > 0, 3C,, > ana 
R = R{t]) > such that the operator from ()4.2p satisfies the following estimate: 

[ l"l<l/9| 

Proof. For |/3| = 0, this follows from [3S1 Lemma 3.3]. We then explain how to 
prove the case with \j3\ > 0. To do this we split K into two parts. We apply p.5p 
and p.Sp with /2 = a/J and fi = hi as in the definition of K from (j4.2|) . etc. 

For the term as in p.8|) . we deduce that it is bounded by a different linear kernel 
with V0j^{q") and Vfi,J^(q). Such a new linear operator has the same property as 
the original K in the center of mass frame with V0j'i{q") and V0j^{q). So the 
Lemma follows from the result for K without derivatives in [25[ Lemma 3.3]. 

For p.Sp . the term with /3i = /3 is bounded by 



sB{p,q,uj) 



R3xS2 p^q" 



dqdu: ji{q) [\{dphi){p')jHq)\ + \jHp'){dphi){q')\} 
sB{p,q,uj) 



'-dqd^ J^q) \idphi)ip)JHq)\. 



This new linear operator has the same property as the original K . So the lemma 
follows in this case from the result for K in [25l Lemma 3.3]. On the other hand, 
for terms with /3i < /?, by compact Sobolev imbedding, it suffices to consider the 



case where \p\ is large, for which the fast decay factor e sIpI^^" in (j3.5p provides 
the small constant ry and we complete the proof. □ 



5. Global solution and rapid decay 

In this final section, we explain how to use the new estimates from the previ- 
ous sections to prove the global existence and rapid decay of nearby Maxwellian 
classical solutions to the relativistic Vlasov-Maxwell-Boltzmann system (|1.6p and 
(|1.9p with (jl.lOp . We prove global existence following the approach from pTlfT^ . 
The decay follows as in the method described in [5D]. Since several of these previ- 
ously elucidated details [TTlIIUlin] are similar, we will simply write down the main 
steps and refer to the prior results for an elaboration of the full argument. We aim 
to make our argument completely precise in the sense that we refer to the exact 
argument which is needed from the previous work [T7l[T9l[20] . 

First in Section lOl we explain the local existence argument. Then in Section [5^ 
we exposit the argument for proving the crucial positivity of the linearized collision 
operator for solutions to the relativistic Vlasov-Maxwell-Boltzmann system (|1.6p 
and (|1.9p with (|1.10p . Finally in Section [5751 we explain how these estimates can be 
used to prove the global in time existence and rapid decay. 

5.1. Local solutions. We now sketch the procedure for obtaining a unique local- 
in time solution to the relativistic Vlasov-Maxwell-Boltzmann system (11.61) . (jl.9p . 
(jl.lOp . These arguments are rather standard [TT l fTQ ] . 

Given a solution [f{t,x,p),E{t,x),B{t,x)] to the relativistic Vlasov-Maxwell- 
Boltzmann system we recall the definition of the instant energy functional £N,i{t) 
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in ()2.ip . We furthermore define the dissipation rate 'DN,i{t) as {(. > 0) 

|7| + I/3|<JV 

We state the following local existence theorem. 

Theorem 3. Fix £>0. 3Mo > 0, T* > such that ifT* < Mo/2 and 

SnAO) < ^^o/2, 

then there is a unique solution [f{t,x,p),E{t,x),B{t,x)] to the relativistic Vlasov- 
Maxwell-Boltzmann system (|L6)) . (fTO)) and (|LT0)) on [0,T*) x Ti^ x such that 

sup <^ SNAt) + ds VnAs) f < Mo- 
0<t<T' I Jo ) 

The high order energy norm E^A^) continuous over [0,T*). 

If Fo{x,p) = J + Ji/^o > 0, then F{t,x,p) = J + .A''^ f{t,x,p) > 0. Further- 
more, the conservation laws il.ll]) . I11.12\) . and I11.13\) hold for all < t < T* if 
they are valid initially (at t ~ 0). 

This local existence theorem can be proven in the standard way using the esti- 
mates in this paper combined with the local existence proof in |19j for the relativistic 
Landau-Maxwell system. The positivity of solutions follows from the proof in [17] . 

5.2. Positivity of L. In this subsection we elucidate the positivity of the linearized 
operator (|1.8|) . L, for any small amplitude solution [f{t, x^p), E(t, x), B{t, x)] to the 
full relativistic Vlasov-Maxwell-Boltzmann system ()1.6p , ()1.9|) and (|1.10p . 
Our main result in this section is as follows. 

Theorem 4. Let [f{t, x,p),E{t, x), B{t, x)] be a classical solution to (1.6\) and lll.9\) 
satisfying M.IO]) . U.ll]) . and \1.13]) . Then there exists an Mq > and a 

(5o = Sn{l\Io) > such that if N > i and 

(5.1) Upvwip + WEm' + WBit)\A < mo, 

h\<N ^ 

then 

Y {Ldy{t),dy{t))>So J2 Wdvm'- 

h\<N h\<N 

In the rest of this section, we always work exclusively with a classical solution 
[f{t,x,p),E{t,x),B{t,x)] to (|1.6p and (|1.9p . This argument proceeds, as is cus- 
tomary, via a careful study of the the six dimensional null space of L, for any fixed 
{t,x); this null space is given by (1 < i < 3) 

(5.2) iV(i) "'=^'span|[\/j,0],[0,VJ],b,Vj,p,Vj], . 
Define the orthogonal projection from L^(Kp) into N{L) by P. Then decompose 

/ = P/ + {I - P}/. 

It is now standard to call P/ = [P_|-/, P_/] e the "hydrodynamic part" of / 
and {I — P}/ = [{I — P}+/, {I — P}-/] the "microscopic part." By separating its 
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linear and nonlinear part, and using L± (P/) = 0, we can express the hydrodynamic 
part of / through the microscopic part up to a higher order term h{f): 

(5.3) + ^ • P±/ T • ^} ^ = '±(^1 - P>^) + '^i^/)- 
This is a decomposition of (|1.6p (with normalized constants), where 

(5.4) ?±({I - P}/) -[dt + ^- V.j {I - P}±/ - i± ({I - P}/) , 

(5.5) /i±(/) T(^E + ^xBy Vp/± ± |i? . iij /± + r±(/, /). 
We further expand P±/ as a linear combination of the basis in (|5.2p as 

(5.6) P±/= |a±(i,a:)+^6,(t,a:)p, +c(t,x)/| 

The positivity of L is obtained via a careful study of the relativistic system of 
macroscopic equations (|5.7p - (|5.1ip : this system was derived in [19]. 

We will sketch the derivation of (|5.7p - (|5.1ip for the convenience of the reader. 
Expand the left side of (|5.3p with respect to the terms in (|5.6p as 

|a°a± + ^ {9^a± tE,} + + P, {d%j + 8^ c) + pOaOcj v/j(^. 

Here 9° = 9* and = d^j ■ For fixed {t, x), this is an expansion of left side of ()5.3p 
with respect to the basis of {efe}, whose components are given by (1 < i,j < 3) 

[Vj, 0], [0, VJ], [pj%/j/p°,0], [0,p,Vj//],p,Vj[1, 1],p,p,\/J//[1, 1],p"\/J[1, 1]. 

We expand the right side of (|5.3p with respect to the same basis and compare the 
coefficients on both sides to obtain the macroscopic equations: 

(5.7) 5"c = h + h,, 

(5.8) d'c + d°h^l, + \, 

(5.9) (1 - 5^j)d'b, + d^h - k, + h,j, 

(5.10) a*a± T-B, = +/i„±, 

(5.11) d''a±=la± + ha±. 

To ease the notation wc define the following index set 



=' {c, i, ij, ai±, a ± | i, j = 1, 2, 3} . 

Thus Ai is the collection of all indices in the macroscopic equations. For X E A4 
each l\{t,x) are the coefficients of ^({I — P}/) with respect to the basis elements 
{cfc}; similarly for each h\(t, x). Precisely, a given l\ can be expressed as 

h=Y. (^'kiihiH - P}/), /-({I - P}/)], efc), e M. 

fe 

Also the h\{t,x) can be computed similarly. 

From (|5.6p and (|1.9p (with normalized constants) wc see that 



(5.12) 



J^^ f [p^/p\ -pVl/p"] ■ {I - P}fdp. 
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Next, we explain how to estimate each of these terms. 
Proposition 9. Suppose that N > 4, then we have 

E ( E ii^'^^ii + ii^'^ii ) E iiii - 

|7|<Ar-l \\eM ) |7|<A' 

This Proposition [9] is proven for the relativistic Landau-Maxwell system in [T9l 
Lemma 11]. Because the structure of the operator l\ is similar, the proof of Propo- 
sition [S] is exactly the same as [1^1 Lemma 11]. The only difference between these 
cases is for the term L± in l\. since the operator L± is in fact very different. However 
because of p.Sp the required estimate for — P}9'''/, Cfc) used in [121 Lemma 
11] is supplied by our Theorem[2l We then refer to [191 Lemma 11] for the rest. 

We now estimate coefficients of the higher order term from (j5.5p . 

Proposition 10. Let \5.1^ he valid for some Mq > 0. Then 

E E \\d'h4<c^o E ii^vii- 

|7|<JVAeA^ |7|<Af 

As in the previous proposition. Proposition [10] is proven for the relativistic 
Landau-Maxwell system in |19[ Lemma 12]. Again the structure of hx is simi- 
lar in both cases, meaning that the proof of Proposition [TUl the same as [THl Lemma 
12]. The only difference between these cases is for the term T± in h^; the operator 
F-t is in again quite different. However the needed estimate for (9T(/, /), e^.) used 
in [ini Lemma 12] is supplied by our Theorem [21 Otherwise the proof is exactly 
the same, and for the full details we refer to the proof in [191 Lemma 12]. 

Next we estimate the electromagnetic field [E{t, x), B{t, x)] in terms of f{t, x,p) 
through the macroscopic equation (|5.10|) and the Maxwell system (|5.12p . 

Proposition 11. Let [f{t,x,p),E{t,x),B{t,x)] be the solution to il.6\} . U.9\) and 
U.10\) constructed in Theorem [3 with the constant (|1.3|) . Let the small amplitude 
assumption i5.1\) be valid for some < Mq < 1 . Then we have 

E {WE{t)\\ + \\d''{B{t)-B}\\)< E wdvm- 

|7|<7V-1 |7|<A' 

Similar to the previous propositions, this Proposition 1 1 II is proven in exactly the 
same way as |19[ Lemma 13] except that we replace the estimates in the proof of 
[191 Lemma 13] with their corresponding analogues herein. Specifically, we follow 
directly the proof of [19| Lemma 13] however we replace the use of [191 Lemma 11 
and Lemma 12] with Propositions 1^ and [TUl respectively. 

Collecting the previous estimates in this section, we can now prove the crucial 
positivity of L from (|1.8p , as stated in Theorem 31 for a small amplitude solution 
[f{t, x,p), E{t, x),B{t, x)]. Once again the proof is the same as the analogous proof 
for the relativistic Landau-Maxwell system from [T^l Theorem 2]. We need only 
replace the estimates used in [191 Theorem 2] with their analogues in this section; 
in particular we replace [T^l Lemma 13] with our Proposition [TTl jl9l Lemma 12] 
with Proposition [TUl and [TUl Lemma 11] with Proposition |U1 We then refer to [TU] 
since otherwise the details are exactly the same. This completes Theorem [H 
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5.3. Global Solutions. In this section we establish Theorem [T] We first explain 
how to derive a refined energy estimate. We use the instant energy functional 

E E Hfmh+ E \\d''m),Bit)r. 

|/3|<,n|7|<Ar-|/J| |7|<Ar 

We also define the refined dissipation rate as 

\f3\<-nih\<N-\l3\ 

Here < m < A^. In these spaces we have the estimate: 

Proposition 12. Fix £ > 0. Let [f{t,x,p),E{t,x),B{t,x)] be the unique solution 
constructed in Theorem\^ which also satisfies the conservation laws il-ll]) . hl.l'^) 
and hl.l3\) . Let the small amplitude assumption liS.l]) be valid. For any given 
< m < N and \I3\ < to, there are constants ^ > and 6m,e > such that 

-J^£m.,eit) + 5m,lDm,l{t) < C*y^ j\J £ N ^{t) 'DN,l{t). 

We point out that the proof of Proposition [12] is exactly the same as the corre- 
sponding proof in [TOl Lemma 14] for the relativistic Landau-Maxwell system. The 
differences are that we use the relativistic Boltzmann estimates from this paper, 
instead of the corresponding estimates from [TH| , and secondly that we include the 
weight £ > 0. For the estimates, we specifically replace [111 Theorem 4] with our 
Theorem [21 [191 Lemma 13] with Proposition [Til [111 Lemma 7] with Propositions 
[71 and [H and otherwise the argument follows exactly the proof of [El Lemma 14] . 
To include the weights ^ > 0, we refer to the argument used to prove [71 Eq (4.6)]. 

Finally we prove the global existence of solutions to the relativistic Vlasov- 
Maxwell-Boltzmann system (|1.6p and (|1.9p . Notice that using the estimates above, 
in particular Proposition [121 this follows from the standard continuity argument as 
in for example [T71ll9|. Thus we have proven all of Theorem [H except for the decay 
rates. But the decay rates in this case follow directly using the interpolation pro- 
cedure from [30] (which was applied to the Newtonian Vlasov-Maxwell-Boltzmann 
system and the relativistic Landau-Maxwell system in [20|). Note precisely that for 
the rapid decay, we use the proof from |20[ Section 2] combined with the differential 
inequality from Proposition [12] in this paper. Q.E.D. 

Acknowledgments. We would like to thank the referees for their careful reading 
of the paper and constructive comments which helped to improve the presentation. 
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